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ABSTRACT 

We present a new, magnetohydrodynamic mechanism for inflation of close-in giant extrasolar plan- 
ets. The idea behind the mechanism is that current, which is induced through interaction of atmo- 
spheric winds and the planetary magnetic field, results in significant Ohmic dissipation of energy in 
the interior. We develop an analytical model for computation of interior Ohmic dissipation, with a 
simplified treatment of the atmosphere. We apply our model to HD209458b, Tres-4b and HD189733b. 
With conservative assumptions for wind speed and field strength, our model predicts a generated 
power that appears to be large enough to maintain the transit radii, opening an unexplored avenue 
towards solving a decade-old puzzle of extrasolar gas giant radius anomalies. 

Subject headings: planets and satellites: interiors — magnetohydrodynamics — methods: analytical 



1. INTRODUCTION 

The detec tion of the first transiting extrasolar planet 
HD209458b ([Charbonneau et al.|2000||Henry et al.|2000 1 
marked the first observation ot a planet whose radius is 



planet 

anomalously large. With the current aggregate of tran- 
siting planets exceeding 60, over-inflated "hot Jupiters" 
are now known to be common (Fig.l), and understand- 
ing their radii has become recogniz ed as an outstandin g 
problem in planetary astrophysics (Baraffe et al. 2010). 



Most proposed explanations require an interior power 
source that would replace the radiated heat from grav- 
itational contraction and cause a planet to reach ther- 
mal equilibrium with a larger-than-expected radius. In 
the context of such solutions, the generated heat must 
be deposited into the interior envelope, i.e. below the 
radiative/convective boundary, in order to maintain the 
core entropy (and theref ore the radius) of the plan et. No- 
tably, eccentricity tides ('Bodenheimer et al. 2001) , obliq- 
uity tides of a Cassini state ( .Winn fc Holman 2005 ) , and 
deposition of kinetic energy to adia batic depths by dy- 
namical and convective instabilities ( Guillot & Showman 



2002 1 have been invoked to provide an extra power source 
in the interior of the planet. It has been shown that th e 
required powers are rather modest ( Burrows et al.|2007 ), 
but it is unlikely that any of the proposed s olutions alone 
are able to account for all observe d radii (Baraffe et al. 
2010| [Fortney fc Nettelmann||2009[ ) 



Here we show that the anomalous sizes of close-in exo- 
planets can be explained by a magnetohydrodynamic 
mechanism. The interactions of zonal winds with the 
expected planetary magnetic field in a thermally ionized 
atmosphere induce an emf that drives electrical currents 
into the interior. These currents dissipate Ohmically and 
thus maintain the interior entropy of the planet. The 
primary controlling factors in our model are the atmo- 
spheric temperature, wind velocity and strength of the 
magnetic field, as they dictate how much current is al- 
lowed to penetrate the interior. Other variables, such 
as metalicity also contribute, but to a smaller degree. 
Our results predict that interior heating of this kind oc- 
curs in all close-in exoplanets with magnetic fields, but it 
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Fig. 1. — Scatter-plot of mass vs. radius of transiting Jovian exo- 
planets. The three planets considered in the text as well as Jupiter 
& Saturn are labeled. The two lines represent the theoretical mass- 
radius relationships for a core-less planet (d ashed) and one with a 
4OM0 core (solid) from |Bodenheimer et ah] ( )2003^ . Planets above 
the dashed line require an mtlation mechanisrn to halt gravitational 
contraction. 

is negligible if the atmospheric temperature is not high 
enough for sufficient thermal ionization to take place. 
Smaller, but hot exoplanets are attributed to heavy el- 
ement enrichment in the interior. While the inflation 
mechanism we present here is general, the quantitative 
modeling in this work is speciflc to HD209458b, Tres-4b, 
and HD189733b which are arguably the better studied 
transiting exoplanets. 

2. STRUCTURAL MODEL & ELECTRICAL 
CONDUCTIVITY 

Unlike Jupiter and Saturn, close-in extra-solar gas gi- 
ants are exposed to high irradiation due to their prox- 
imity to parent stars. This forces their atmospheric 
temperature-pressure profiles to be signific antly shal- 
lower than their solar system counterparts (Fortney & 
Nettelmann 2009). In particular, the lower atmospheres 
{F > 0.1 Bars) of hot Jupiters are believed to be almost 
isothermal while the radiative/convective boundaries are 
though t to lie at P ~ ICIO — 10 00 Bars, depending on the 
planet ( Showman et al.||2008| ). 

The isothermal sections or extrasolar gas gian t atmo- 
spheres often reach temperatures close to 2000 K ( Spiegel 
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et al. 1 2009 ) and in some cases, even higher (Borucki et al. 
2009| . 'iTiese temperatures are not high enough to ionize 
H or He significantly, however, alkah metals such as Na 
and K will be partially ionized. As a result, electrical 
conductivity in the interior of a hot Jupiter is dominated 
by ionization of hydrogen, while in the outer region of 
the planet, electrical conductivity is primarily due to the 
ionization of alkali metals, with the transition between 
the two inoization regimes taking place at P ~ 300 Bars. 
Thermal ionization is governed by the Saha equation: 



ntrie 



exp (-Ij/kbT) 



(1) 



where nj and rij' are the total and positively ionized 
number densities of constituent j respectively, = 
^ n'j is the total electron number density, me is the elec- 
tron mass, kb is Boltzmann's constant, T is temperature, 
h is Plank's constant, and Ij is the ionization potential 
of constituent j. If the ionization is far from complete 
(n^ ^ fju), the abundances of alkali metals, fj are held 
constant, and the atmosphere is isothermal, it is easy to 
show that the electron number density takes on an ex- 
ponential profile with an ionization scale-height that is 
twice as large as the density scale-heght: 



no 
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where x's are the RHS's of equation (1), vq is the ra- 
dial distance at some reference point (P = 10 Bars) and 
H ~ kbT/ is the density scale-height. In our ionization 
calculations, we considered the following alkali metals: 
Na, K, Li, Rb, Fe, Cs and Ca. Their abundances and 



ioni zation potentials were inferred from Lodders ( 1999 1 
and Cox et al. ( 2000[ ) respectively. 

The atmospheric temperatures above the isothermal 
layer differ significantly from planet to planet. In partic- 
ular, thermal inversions have been detected in the atmo- 



spheres of HD209458b (Burrows et al..2007J and Tres-4b 



(Knutson et al. 2009) but not in HJJ189733b. In our 



models, we ad opt atmospher ic temperature profiles sim- 
ilar to that of 
Tres-4b, and t 
HD189733b. 



Spiegel et al. (2009) for HD209458b and 



re ID profile of |Fortney et al. (2010) for 
The relatively cool temperatures attained 
above P < 0.1 Bars are of significant importance to our 
models because they provide insulating shells which are 
impenetrable to radial current. Consequently, current 
loops are necessarily setup through the interior, and any 
current fiowing in the ionosphere is not relevant. We 
place the radiative/convective boundary at P ^ 100 Bars 
in all of our models. 

We did not have to explicitly compute the ionization 
fractions of H and He, as they are published in the equa- 
tion of state ( Saumon et al.|1995 ), which we employed in 
our model. In particular, we used the "interpolated" 
version of the equation of state, where ionization oc- 
curs smoothly with pressure and temperature. Although 
the planetary structure was core-less, we mimicked the 
presence of a core by c hanging the Helium c ontent from 
Y = 0.24 to F = 0.3 ( [Burrows et~al]|2003[ ) in some of 
our models. 
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Fig. 2. — Electrical conductivity profile of the nominal 
HD209458b model with Tiso = 1700K, Y = 0.24 and Z = Ix solar. 
The inset is a magnification of the profile in the outer part of planet. 
The model parameters R, 5 and 7 are labeled. The dashed lines are 
functional approximations to the conductivity profile where zonal 
fiow is present. The highlighted region corresponds to the upper 
convective envelope (100 — 3000)Bars, where most of the interior 
dissipation takes place. 

Having computed the electron number de nsity, the 
electrical conductivity of a gas is given by (Tipler & 
Llewellyn|[2002| ) 



n m„Ay 8khT 



TTirie 



(3) 



where n and A are the number density, and the number 
density weighted cross-section of everything other than 
electrons. Strictly speaking, the above equation is only 
valid for non-degenerate gas. However, by the point mat- 
ter becomes degenerate in our models, the resistivity is 
completely negligible and the details of its profile have 
no noticeable effect on the results. 

Since we are only interested in the part of the planet, 
interior to the atmospheric temperature minimum, we 
define the model radius r = P as the point of maximal 
conductivity in the atmosphere (P — 75 mbars) , and we 
set the outer edge of our model at the conductivity min- 
imum, r = R + ^ (P = 30 mbars) . We place the bottom 
boundary of the "weather" layer of the atmosphere at a 
pressure of P = 10 Bars and denote it as r = R — 6. Con- 
sequently, the "inert" layer of the atmosphere is between 
100 ^ P ^ 10 Bars. A computed electrical conductiv- 
ity profile for HD209458b is presented in figure (2), along 
with a simplified conductivity profile resulting from equa- 
tion (2). Because the functional profiles (dashed curve) 
are in good agreement with the numerically computed 
profile, we utilize them in all future calculations (see ap- 
pendix). 

3. ANALYTICAL THEORY 
Global circulation models 



Langton 
have 



& Laughlin 
shown 



that 



Showman et al.'2008; '2009} 
Hau scher, 2009] ) 
specifically 



2008| [Menou fc 
winds on hot Jupiters, 
HD209458b and HD189733b, can attain velocities of or 
der V ^ Ikm/s. It appears that two qualitative wind 
patterns are present. In the upper atmosphere (P < 
30mbars), wind fiows from the sub-stellar point to the 
anti-stellar point symmetrically across the terminator. 
Deeper down, a strong eastward zonal jet develops. Im- 
portantly, the development of zonal .je ts have been ob 
served in virtually all simulations (see 
(j2009j) for a comprehensive review). 



Showman et al. 
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Prior to obtaining a formal solution to the problem, 
we can identify some of its features. First, if the planet's 
dipole moment is aligned with the rotation axis and we 
consider only zonal flow, then there is azimuthal symme- 
try. Second, it is immediately apparent from the geome- 
try of the zonal jet and the dipole field that the induced 
current will be meridional. In the atmosphere, we expect 
the current to flow from the poles to the equator where 
it penetrates the interior of the planet and completes the 
loop (Fig 3). 

The general induction equation can be written as: 



dB 



VxAVxB +Vx Ivx B 



(4) 



where B is the magnetic field a nd A = l//iocr is the mag- 
netic diffusivity (Moffatt 1978 1 . We express the magnetic 
field as a dipole background component and an induced 

component: B = B dip + Bind with V x Bdip — 0. This as- 
sumes no dynamo generation in the region. The induced 
magnetic field will tend to point in the same direction 
as the velocity field, so we can make the approximation 

V X B m V X Bdip- We assume that the prescribed veloc- 
ity field and background magnetic field are not strongly 
modified by the induced field i.e. Rm = vL/X < 1, an as- 
sumption satisfied in our models with T^ 1700K. Finally, 
we seek a steady-state solution, so we require dB/dt — 0. 
With these assumptions, the induction equation simpli- 
fies to: 



V X A V X B. 



= V X 



(v X Bdip^ ■ 



(5) 



We can "uncurl" this equation and use Ampere's law 
V X B = hqJ to recover Ohm's law: 



(6) 



By continuity, V • J must vanish. As a result, 

^ • f7^$ -adfx '^drp) ■ (7) 



If the conductivity takes on an exponential form, there 
exists an analytical solution for $ and in our models, we 
confine the atmospheric fiow to the region where conduc- 
tivity is exponential. In the interior region, the electric 
potential is also governed by the above equation, with the 
right-hand side is set to zero. However, since the interior 
conductivity does not take on a simple analytical form, 
the above equation there must be solved numerically. 

We take a nominal value for the "strength" of the field 
at the surface of the planets to be \\B\\ji — 10~'^T, ap- 
proximately the value expected from scaling the field via 
the Elsasser number A = aB^ /2pil ^ 1, where 57 is the 
planetary rotation rate (assumed tidally locked). The 
magnetic field scaling argum ent based on energy flu x 



also suggests a similar value (Christensen et al. 
For comparison, Jupiter's surface field is \\B\ 



2009) 



4.2x10 ''T ( Stevenson 2003 1 . We approximate the zonal 
wind as V (X Vm sai{6)(j) where Vm is the maximum speed 
attained by the wind and set Vm — 1 km/s (see appendix 
for more details). 

Once we have the solution for the current, we can com- 
pute the total Ohmic dissipation rate below some radius 




Fig. 3. — Side view cross-section of induced current due to zonal 
wind flow. The interior vector field, plotted with small arrows, 
is a quantitative result of the model. The large semi-transparent 
arrows are illustrations. The yellow shell in the inset represent 
the region to which we confine the zonal flow (10-0.03 Bars). The 
orange region denotes the region of interior heating. 



cr(r) 



dV. 



(8) 



In order to satisfy continuity, the magnitude of the cur- 
rent density must be constant along its path in the inte- 
rior. As a result, it is apparent from the above equation 
that most of the dissipation takes place in the upper lay- 
ers of the planet, where conductivity is not too great, and 
the solution is insensitive to the details of the conductiv- 
ity profile in the deep interior, as long as it remains high. 
The Ohmic heat that is generated in the convective enve- 
lope of the planet replaces gravitational contraction, and 
is lost by radiative cooling at the radiative/convective 
boundary. Consequently, to ensure a null secular cool- 
ing rate, we need the Ohmic dissipation rate to at least 
compensate for the the ra diative heat flu x at the radia- 
tive/convective boundary (Clayton|l968|. 



4. MODEL RESULTS 

It has been shown that extra-solar gas giants require 
between 10~^ and 10~^ of the irradiation they receive to 
be de posited into the adiabatic interior to maintain their 



radii ( Bodenheimer et al. 2001 Burrows et al. 2007 Ibgui 



et al. 



2U09[ ), although the exact number depends on the 



metaliicity of the atmosphere and the mass of the heavy 
element core in the interior of the planet. Under the as- 
sumption of solar metaliicity and no core, IID209458b 
requires 3.9 x lO^^W, Tres-4b requires 8. 06 x 10^°W and 
HD189733b requires no heating at all (Burrows et al 
[2007; Ibgui et al. 2009). Within the context ot our model 
HD209458b and HD189733b are easily explained. To 
adequately explain Tres-4b however, we require an en- 
hanced (lOxsolar) metaliicity in the atmosphere to re- 
duce the required heating down to 5.37 x lO^^W. 

Table (1) presents a series of models with various tem- 
peratures. Helium contents, and metallicities of the plan- 
ets under consideration. Upon inspection, it is apparent 
that the global heating rate scales exponentially with 
temperature, and as a square root of the metaliicity. 
Both of these scalings can be easily understood by noting 
that scaling the conductivity profile by a multiplicative 
factor causes a corresponding change in dissipation while 
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equations (1) and (2) relate temperature and metallicity 
(i.e. /) to the conductivity. 

It is also noteworthy that the models with a simulated 
core produce approximately the same amount of heating 
as the coreless models. This is because most of the dis- 
sipation takes place in a region where ionization of alkali 
metals still dominates the electrical conductivity and the 
somewhat hotter interior isentrope makes little difference 
- recall that the value of the conductivity is unimportant 
in the deep interior. 

There are a number of other scalings present within our 
model. For instance, the total Ohmic dissipation rate is 
proportional to the squares of the wind speed and the 
strength of the magnetic field, P cx (i?/10^T)^(?;m/lkm 
s~^)^. Additionally, to leading order, the dissipation in 
the atmosphere scales linearly with the thickness of the 
atmosphere, while the interior dissipation approximately 
scales quadratically. Consequently, along with the con- 
ductivity effect, hotter atmospheres also lead to more 
dissipation by virtue of a physically larger atmosphere. 

It is important to understand that Ohmic heating does 
not only affect the interior. Because the induced cur- 
rent's "return path" lies in the atmosphere (Fig. 3), the 
atmosphere also gets heated. This heating, along with 
magnetic drag on the flow are the limiting factors of our 
theory. 

Consider the nominal case of HD209458b with Tiso = 
1700K and Z = Ix solar. In this model, the heating re- 
quired to inflate the planet is attained at a depth of ~ 90 
Bars, essentially at the radiative/convective boundary. 
Ohmic heating in the atmosphere is small (only ~ 4%) 
in comparison with the irradiation, having little effect 
on the planet's evolution or structure and i?m ~ 0.3. 
In other words, the assumptions implicit in our calcu- 
lation are satisfied and the mechanism seems to explain 
the transit radius adequately. However, if we go to the 
model with Tiso — 2000K, ~ 3, the Ohmic dissipa- 
tion in the atmosphere is comparable with the insolation, 
and the assumptions of our model may no longer apply. 

The nominal model of Tres-4b with T^so = 2250K also 
runs into the same problem. Here, Rm ~ 15, and the 
Ohmic dissipation in the atmosphere is again comparable 
with the insolation. However, if we imagine that mag- 
netic drag reduces the wind velocity by a factor of ~ 3, 
our results fall in the right ballpark to explain Tres-4b's 
radius, in the scenario where its atmospheric opacity is 
super-solar. Finally, consider the model of HD189733b. 
For this configuration of parameters, our mechanism does 
not predict a significant amount of Ohmic dissipation 
at adiabatic depths, consistent with an un-inflated ra- 
dius. A similar scenario is observed for the model of 
HD209458b with Ti,o = 1400K. Overall, it appears that 
within the current setup of the model, the cumulative 
heating below the weather layer of the atmosphere i.e. 



r < R— S is of order a fewx lO^'^ of the heating that takes 
place in the atmosphere. Provided that this ratio of mag- 
nitudes holds up in a more dynamical treatment of the 
problem, it can provide an upper limit to the maximum 
inflation that can be explained with Ohmic dissipation. 



5. DISCUSSION 

In this letter, we have presented a new, magnetohy- 
drodynamic mechanism for inflation of extrasolar gas gi- 
ants. Our calculations show that the heating, necessary 
to maintain the seemingly anomalous radii of transiting 
exo-planets, naturally emerges from considerations of in- 
teractions between partially ionized winds and the plane- 
tary magnetic field. Interestingly, there seems to be a set 
limit to the extent that Ohmic dissipation can heat the 
interior, making this theory testable. Currently, there 
is significant uncertainty with respect to the calculation 
of the required interior heating, because core masses are 
unknown. Ho wever, dynamical det e rminations of inte- 
rior s tructure (Batygin et al. 2009 Ragozzine & Wolf 
2009 1 may allow us to resolve the degeneracy for a frac 



tion of observed planets, and provide a solid test-bed for 
the mechanism we've presented here. 

There is a number of interesting additional questions 
that our model inevitably brings up. First, recall that 
our treatment of the induction equation is kinematic. In 
reality, flow modification by the Lorenz force may play 
an important role in determining the actual wind pat- 
terns. While this effect may be small for HD209458b 
and HD189733b, weather on hotter planets, such as Tres- 
4b or Wasp-12b may be more intimately linked with 
their magnetic fields, calling for a magentohydrodynamic 
treatment of the atmospheric circulation. Generally, 
when zonal winds interact with a background dipole field, 
they give rise to poloidal current which in turn, gives rise 
to a predominantly toroidal, unobservable field. How- 
ever, the dayside-to-nightside flows that are present at 
higher levels in the atmosphere may modify the flow in 
an interesting way that may eventually be astronomically 
observable. 

Second, we are neglecting the stellar magnetic field. 
The star's magnetic field is likely to be considerably 
smaller than the planetary field at the planetary orbital 
radius, but induction by stellar field as well as coupling 
of the stellar and planetary magnetic field lines is cer- 
tainly plausible. This too, may produce an astronomi- 
cally observable signature. Finally, we are neglecting the 
effects the induced current in the interior on the plan- 
etary dynamo. Considerations of this sort may influ- 
ence the background magnetic field of the planet. All of 
these aspects call for a self-consistent treatment of the 
full problem. Such calculations would no-doubt provide 
further insight into the physical structure of extra-solar 



APPENDIX 

We approximate the electrical conductivity profile in the atmosphere with exponential functions: 



age 



R-S <ri^ R 
R<r s^R + -f 



(1) 



where ag and a-y are the conductivities at r = R — S and r = R respectively, while Hg and Hj are the conductivity 
scale-heights in the corresponding regions. We prescribe a parabolic radial dependence to the zonal fiow over the 
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TABLE 1 

OhMIC dissipation acquired at various pressures in various models op HD209458B, TRES-4B, AND HD189733B 



Planet 


Y 


Ttso (K) 


Z (x solar) 


P |P < 10 BarsJ (W) 


P [P > 10 BarsJ (W) 


P > 100 BarsJ (W) 


HD209458b 


0.24 


1400 


1 


2.30 X 


lo^y 


2.23 X 


lOi^ 


1.09 X 


10i« 


HD209458b 


0.24 


1400 


10 


7.28 X 


1019 


7.06 X 


10" 


3.43 X 


1016 


HD209458b 


0.24 


1700 


1 


1.14 X 


1021 


1.01 X 


1019 


5.60 X 


101^7 


HD209458b 


0.24 


1700 


10 


3.61 X 


1021 


3.19 X 


10i9 


1.77 X 


IQi* 


HD209458b 


0.24 


2000 


1 


1.22 X 


1022 


3.24 X 


1020 


7.09 X 


1019 


HD209458b 


0.24 


2000 


10 


3.89 X 


1022 


1.05 X 


102i 


2.29 X 


102o 


HD209458b 


0.3 


1400 


1 


2.22 X 


1019 


1.30 X 


lOi^ 


9.18 X 


1014 


HD209458b 


0.3 


1400 


10 


7.01 X 


1019 


4.10 X 


1017 


2.89 X 


IQis 


HD209458b 


0.3 


1700 


1 


6.97 X 


1020 


7.67 X 


lOis 


8.02 X 


1017 


HD209458b 


0.3 


1700 


10 


2.21 X 


1021 


2.43 X 


1019 


1.90 X 


10i« 


HD209458b 


0.3 


2000 


1 


1.38 X 


1022 


3.13 X 


1020 


4.05 X 


1019 


HD209458b 


0.3 


2000 


10 


4.52 X 


1022 


1.05 X 


1021 


9.42 X 


1019 


Tres-4b 


0.24 


2000 


1 


6.87 X 


1022 


2.57 X 


1021 


1.42 X 


1020 


Trcs-4b 


0.24 


2250 


1 


1.44 X 


Iq23 


3.33 X 


1021 


3.68 X 


1020 


Trcs-4b 


0.24 


2500 


1 


4.62 X 


Iq23 


7.87 X 


1021 


1.54 X 


1021 


Trcs-4b 


0.3 


2000 


1 


4.80 X 


1022 


9.56 X 


1020 


3.16 X 


1019 


Trcs-4b 


0.3 


2250 


1 


1.98 X 


Iq23 


5.92 X 


1021 


6.16 X 


1020 


Trcs-4b 


0.3 


2500 


1 


5.13 X 


1023 


8.75 X 


1021 


1.55 X 


1021 


HD189733b 


0.3 


1500 


1 


9.94 X 


IQlS 


2.65 X 


1016 


1.00 X 


10" 



thickness of the atmosphere, 6, and maintain the velocity constant over the outermost thin shell, 7: 

0<r^R-S 
^Vmsm9$ R-SKr^R + j 

where 

[1 R<r^R + ^ 

Assuming alignment of the dipole moment and the rotation axis, the background dipole magnetic field can be expressed 
as follows: 

Bdip = fxk (^-^) I (4) 



With these expressions, we can decompose the angular part oivxB into spherical harmonics. Upon inspection, one 
finds that the only harmonic of interest has £ = 2 and m = 0. Consequently, we write the potential as $ = g{r)Y2{9, cj)) 
and equation (7) becomes a scalar equation. 

Because the outer edge of our models is set at an insulating shell, we require the radial current at r = i? + 7 to be 
zero: 

5.(i? + 7) = ^5^(^^. (5) 

This boundary condition is appropriate when the electrical resistance that the current will encounter radially, greatly 

exceeds that of a path confined to a surface i.e. fl^^'^'^ a~^dr ^ -R/g^ a^^dO. This criterion is satisfied in our models. 

With this boundary condition, the radial part of the solution to equation (7) in the outermost shell {R < r ^ R + "f) 
reads: 

R+r+~t 

e 

" 90if^r3 (1 - AH^ + m^) {12H^ - m^{R + 7) + (i? + 7)2) 

X {l2^/^kvme^^'^ (1 - AH^ + 6iJ^) {6H'^{2R - 5r + 27) 

- 2H^{R'^ - 3Rr - 3r^ + 2i?7 - Srj + 7^) - rCr^ + {R + jf)) 

- SOe^H^iAH^ + r){12H^ - 6H^{R + 7) + (i? + 7)2)^1 

+ 5e^ {2AH^ - ISH^r + 6H^r^ - r^) {l2H^ + 6Hj{R + j) + {R + 7)^) ^1), (6) 

where Ai is an undetermined constant of integration. In a similar fashion, we can write down the solution to the radial 
part of equation (7) in the region {R — S < r ^ R): 

1 5^2 {-24HI + 18Hlr - GHgr^ + r^) bAsHge' * (AHs + r) 
9s{r) = T^[ ^TEH 



15r3 ' -AHg + l - AHg + 1 
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12//|(4iJ5 + r)Ei ( - 



X (-192i/| + '^mlr + 2(-24ij| + 18i/|r - ei/^-r^ + r^) log(r) 

- 2Hsil2r^ + - 5^) + r{4rR - + (5^))), (7) 

where A2 and A3 are again undetermined constants, and Ei is an exponential integral: Ei(x) — T*^^- Although 
equation (7) must be solved numerically in the interior, towards the center of the planet, where conductivity can 
be taken to be constant, it reduces to Laplace's equation. As a result we can use the polynomial eigenfunction 
ginti'r) = A^r"^ in the vicinity of the origin, and is the last undetermined constant. The four constants of integration 
are determined by continuity. 
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